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1 Introduction 



In this paper we continue our discussion on the Coulomb branches of gauge theories, which admit an 
exact treatment in the low-energy sector |l| and are related to integrable systems [||-[[l5|. These are 
the theories with eight real supercharges (in accordance with common terminology that means N = 2 
SUSY in four dimensions and Af = 1 SUSY in D = 5, 6) with only four noncompact directions to have 
non-trivial instantonic contribution and, therefore, a non-trivial prepotential. It means, in particular, 
that we consider 5d theories with the target space of topology i? 4 x S 1 and 6d theories with the target 
space of topology R 4 x T 2 . Below we give a detailed analysis of the results for the prepotentials of 
these higher dimensional theories. This paper continues our investigation of the prepotentials started 
in |lC| ] and, in particular, gives some support to the integrable structures of higher dimensional theories 
proposed in [|4| [Dl . 

We consider here the SUSY theories with massive hypermultiplets in fundamental representations. 
Therefore, in accordance with || 0, |lq], the corresponding integrable structures are given by (classical) 
spin chains: the XX Z chain in 5d and XYZ (Sklyanin) chain in 6d. Below we concentrate on the issues 
related to the prepotentials and associated analytic structures on Riemann surfaces: differentials, their 
periods etc. In particular, we calculate the perturbative prepotentials which coincide with the expected 
from general considerations of quantum field theory. It gives some support to the identification of the 
integrable spin-chains with the SUSY QCD in higher dimensions. 

An integrable system, i.e. a spectral curve and a generating 1-differential dS Jl^j , allows to define 
the prepotential in the framework of the Seiberg-Witten (SW) anzatz where it can be identified 
with the (logarithm of) tau- function of the Whitham hierarchy In sect. 2 we remind this procedure 
for the standard case of Ad M = 2 QCD. In sect. 3 we extend it to the 5d case, in particular, we 
advocate there the most effective method of calculating the prepotentials based on residue formulas 
first studied in the framework of the SW theory in connection with the associativity equations || [l0| . 
The procedure becomes especially explicit when computing the perturbative prepotentials. In sect. 4 
we discuss 6d theories described by the XYZ spin chain and sect. 5 contains some discussion. Appendix 
contains some definitions and formulas on elliptic functions we use throughout sect. 4. 



2 Ad supersymmetric theories: the Seiberg-Witten anzatz 

First, as a warm-up example we consider the integrable system corresponding to the 4g? SUSY QCD 
and explain how the data obtained from integrable theory can be associated with the objects in 
nonperturbative SUSY theories, in the context of the Seiberg-Witten anzatz jjj. Our presentation 
follows the line worked out in j|, B, ||, [n], [lB) and some details omitted below can be possibly found 
there. 

2.1 Ad theories and XXX chain 

Let us remind, first, the definition of the sl(2) (integrable) XXX spin chain. The most convenient 
formulation is based on the 2x2 Lax matrix (see for example [0), which for the sl(2) XXX chain 
is 

3 

L{\) = \-l + Y,S a -v a . (2.1) 

a=l 



and the Poisson brackets on the space of the dynamical variables S a , a — 1, 2, 3 in the 2x2 formalism 
are implied by quadratic r-matrix relations 



{L(A) ? L(X')} = [r(X - A'), L(A) ® L(A')] , (2.2) 

with the rational r-matrix 

1 3 

r(A) = -^ ff "^" (2.3) 



o=l 



Following standard procedure [EJ, one should consider the c/iam of the Lax operators (2.1) associated 



to each site of the chain and commuting with each other at different sites. The linear problem in the 
spin chain Jl7| can be defined by 

Li(A)*i(A) = * i+ i(A) (2.4) 
where ^(A) is the two-component Baker- Akhiezer function. One can introduce the transfer matrix 



acting from i-th site to i + n-th taking product of the Lax matrices (2.5). The crucial property of the 



relation (2.2) is that it is multiplicative, i.e. any transfer matrix, being product of the form 

T n (X) = L n (X) . . . L\[X) 

(2.5) 

T n (A)*i(A) = * l+ „(A) 

satisfies the same relation 

{T n (A) f r„(A')} = [r(A - A'), T„(A) ® T n (A')] , (2.6) 

provided all L,: in product ( [2~5| ) are independent, {Li,Lj} — for i ^ j. The relations ( pT^ ) are also 
true for the inhomogeneous T-matrix (A — > A + Ai at i-th site; in the applications to M = 2 SUSY 
gauge theories the chain has to be taken of the length N c |j| ) 

T Nc W= II L *(^-\) (2.7) 

Af c >i>l 



still satisfying (|2.6|) with the same r-matrix (2.3) 



In the sZ(2) XXX case the r-matrix relations fl2.2| ) are equivalent just to the common si (2) 
commutation relation (written in terms of the Poisson brackets) 

{S a ,S b } = ie abc S c , (2.8) 

i.e. {S a } plays the role of angular momentum ("classical spin" ) variables giving the name "spin-chains" 
to the whole class of systems. The algebra ( |2.8| ) has an obvious Casimir operator (an invariant, which 
Poisson commutes with all the generators S a ), 

3 
a=l 

so that 

det 2x2 L(A) = X 2 -K 2 , 



N c N c 

det 2x2 T Nc (X) = n det 2x 2 ^(A Xi) = [] ((A - A 4 ) 2 - K 2 

i=l i=l 



(2.10) 



H(X-mt)(X-mr) = Q 2Nc (X), 



i=i 



where we assumed that the values of spin K can be different at different nodes of the chain, and [] 

mf = A l ± Ki. (2.11) 



While the determinant of monodromy matrix (2. 1C ) depends on dynamical variables only through the 
Casimir functions Ki of the Poisson algebra, the dependence of 



Pn c W = ^Tr 2 x 2 7k(A) = J] (A A,) 



(2.12) 



upon dynamical variables is less trivial. Still, as usual for integrable systems 



{P Nc (X),P Nc (X')} = ^{TyT Nc (X),TyT Nc (X / )} = ^Tt [r(X - A'), T Nc (X) ® T Nc (X')] = (2.13) 

i.e. the polynomial P/v c (A) ( |2.12 ) depends on Sa^ only through the Hamiltonians of spin chain 
(Poisson-commuting with each other) . 

The (quasi)periodic boundary conditions imposed to the Baker-Akhiezer function 



*i+iV c (A) = w*i(A) 
give rise to the spectral curve equation for the periodic spin chain 

det (Tjv c (A) - w ) = 0, 

In the particular case of sl(2) XXX spin chains, the spectral equation acquires the form ^ 

QnAX) 



w + 



w 



= 2P Nc (X) 



(2.14) 



(2.15) 



(2.16) 



where the number of matter multiplets Nf < 2N C depending on a particular "degeneracy" of the full 
chain and, as follows from (2.1), ( 2.12 ) A; =0. 

It is relatively easy to get the exact expressions for the Hamiltonians (the explicit examples of 
monodromy matrices and Hamiltonians for N c = 2 and 3 can be found in The Hamiltonians 

depend non-trivially on the A^-parameters (inhomogeneities of the chain) and the coefficients in the 
spectral equation (2.15) depend only on the Hamiltonians and symmetric functions of the mass- 
parameters (2.11), i.e. the dependence of {A^} and {Ki} is rather special. This property is crucial 
for identification of the m-parameters with the masses of the matter supermultiplets in the J\[ = 2 
SQCD. 



2.2 Spectral curves and prepotentials 

Let us turn now to the 4c? SW anzatz [|l] , which can be formulated in the following way. The J\f = 2 
SUSY vector multiplet has necessarily (complex) scalars with the potential V(<f>) = Tr[e/>, <fr\ 2 whose 
minima (after factorization over the gauge group) correspond to the diagonal ([0,0^] = 0), constant 
and (in the theory with SU (N c ) gauge group) traceless matrices. Their invariants 

N c 

det(X-<j>) = P Nc (X)=Y,SN c -kX k (2.17) 

fc=0 

1 Eq. (2.11) implies that the limit of vanishing masses, all = 0, is associated with the homogeneous chain (all 



Xi = 0) and vanishing spins at each site (all Ki = 0). 

2 In general case of sl(p) chains the spectral equation is of the p-th order in w, see |5l La] and below. 



(the total number of algebraically independent ones is rank SU(N C ) = N c — l) parametrize the moduli 
space of the (Coulomb branch of the) theory. Due to the Higgs effect the off-diagonal part of the 
gauge field A M becomes massive, since 



fo,A,Jy = (&-&)Ajf (2.18) 



while the diagonal part, as it follows from (2.18) remains massless, i.e. the gauge group G — SU{N C ) 
breaks down to f7(l) rankG = U{\) N °~ l . 

The effective abelian U(l) Nc ~ 1 theory is formulated in terms of a (finite-dimensional) integrable 
system: the spectral curve £ defined over the genus-dimensional subspace of the full moduli space, 
e.g. 

V(X,w) = 2P N (A) -w- — = (2.19) 

to 

for the pure SU(N C ) gauge theory and the generating differential 

dS = A— (2.20) 

w 

whose basic property is that its derivatives over N c — 1 moduli give rise to holomorphic differentials. 
From the point of view of integrable system, the spectral curve express the condition of the common 
spectrum of two operators, usually the Lax operator C and some auxiliary operator, which was taken 



above to be the shift operator T/v c (2.5). From practical point of view one can work in the basis of 
eigenvalues of one of the operators: in the formalism presented above we have chosen the basis of 
the Lax eigenfunctions - the wave functions which explicitly depend on the Lax eigenvalue A, then 
another operator acts in this basis as a A-dependent matrix. As a result, the full spectral curve £ 
arises as a ramified covering over some bare spectral curve: 

V(X;z) = (2.21) 

where A or z (z = \ogw in ( 2.19Q ) can be considered as a co-ordinate on bare spectral curve So- ln 



the case of the gauge group G — SU{N C ), the function V, caused by (2.17), is a polynomial of degree 



N c in A. If one adds the fundamental matter hypermultiplets [|19| one should also complete the set of 



moduli by their masses m a (a = 1, . . . , Nf) and modify the form of the curve (2.19) to (2.16). Then, 
natural generating differential is 

dS = \— (2.22) 

w 

There is, however, another generating differential usually choosen in literature jl], [jj, |l9), which is 
expressed through the variable W = 7 = . • 

dS = A^ (2.23) 

This variable is naturally associated with the curve written in some more symmetric formal 

W+±= 2PnM (2 24) 

W V /Q^(A) (2 ' 24j 

The both choices of the differentials lead to the same results Jul. In this paper, we use the differential 



(2.24) as better suited for perturbative calculations. 



3 In terms of brane picture it depends on the exact form of embedding of a spectral curve into (eleven-dimensional) 
target space (see, for example, |l2[ |l4[ jf^ , |2c| for more detailed discussion of this point). 



The coefficients si of the function V parametrize the (subspace of) moduli space of complex 
structures Mt. of the curve £. The Hamiltonians or action variables (integrals of motion) are specific 
co-ordinates on moduli space. From the four-dimensional point of view, the co-ordinates sj include 
- (the Schur polynomials of) the adjoint-scalar expectation values hk — j(Ti<p k ) of the vector Af = 2 
supermultiplet, as well - the masses of the hypermultiplets. One associates with each 

handle of S a gauge modulus and with each puncture - a massive hypermultiplet with mass given by 
residue of dS in the puncture. 

The generating 1-form dS = Xdz is meromorphic on E (hereafter the equality modulo total deriva- 
tives is denoted by "="). The prepotential is defined in terms of the cohomological class of dS: 



ai= <b dS, |^ = / dS 

'At da i JBj 



(2.25) 



Ai o Bj = 5 



i.j- 



The cycles Aj include the A^s wrapping around the handles of £ and A a 's, going around the sin- 
gularities of dS. The conjugate contours Bj include the cycles Bi and the non-closed contours B ai 
ending at the singularities of dS (see [[lOj for details). The integrals J g dS are actually divergent, but 
the coefficient of divergent part is equal to residue of dS at particular singularity, i.e. to a a . Thus, 
the divergent contribution to the prepotential is quadratic in a a , while the prepotential is normally 
defined modulo quadratic combination of its arguments (which just fixes the bare coupling constant). 
In particular models § A B dS for some pairs of conjugate contours are identically zero on entire A4: 
such pairs are not included into the set of indices {/}. 

Note that the data (E, dS) with such properties are exactly the definition of the integrable system 
in the sense of jl6| (see also and references therein for details). The period matrix of E: Tij(ai) = 
d ®.g a . as a function of the action variables ai gives the set of coupling constants in the effective abelian 
{/(l)^c- 1 theory while action variables themselves are identified with the masses of the BPS states 
M 2 ~ |na+ ma/)| 2 with the (n, m) "electric" and "magnetic" charges. As an output of all these 
definitions, one can derive the following residue formula J9|, |l^| for the derivatives of prepotential 

d^T dTu duidujjdwK /„ 

= res — (2.2b) 



dajdajdaK dax dz=o dzdX 

where duji are the canonical holomorphic I-differentials satisfying f A dcuj = Sij. Their basic property 
§ B dujj = Tij can be obtained by differentiating the meromorphic differential dS w.r.t. to the moduli 
a/. 

This formula can be effectively used for the calculations of the prepotentials term by term, in 
particular, for the perturbative calculations. 

2.3 Perturbative prepotential 

Now, following mn we discuss the simplest case of degenerate curves corresponding to the perturbative 



formulas. In the perturbative limit ( — > 0) the second term in (2.24) vanishes, the curve acquires 



the form 



(a rational curve with punctures Q) and the generating differential (2.23) turns into 

We remind, first, the simplest case of the pure gauge SU(N C ) Yang-Mills theory ( [2.19| ), ( [2.20 ), i.e. 
Q(A) = 1. Now the set of the N c — 1 independent canonical holomorphic differentials is 



1 1 \ „ _ X iNc dX 

X-Xi X-XnJ ~ (A-A i )(A-A JV J' 



duji = ( r - — ) dX = ^ — , i = 1, N c - I, Xij = Xi - X, (2.29) 



and the ^4-periods at — Xi (the independent ones are, say, with i = 1, . . . , N c — 1) coincide with the 
roots of polynomial P/v c (A) ( [2.12 ). The prepotential can be now computed via residue formulas [jl0| 

dwidajjdwk , . 

Tij k = res - — - J (2.30) 
diogp„ c =o d\ogP Na {X)dX 

which can be explicitly calculated. The only technical trick is that it is easier to compute the residues 
at the poles of (iw's instead of zeroes of dlogP. This can be done immediately since there are no 
contributions from the infinity A = oo. The final results have the following form 



1 6^1 

^ A lfc A Wc ^ A fcJVc 

^ X; AT X^AT ' J XhAT X; AT A™' (2.31) 



i#JV e l^ti,],k,N a 



giving rise to the prepotential formula 



^ = iE/ (4) (a.) 



4 ij (2.32) 
f^\x) =x 2 \ogx 2 

Now let us turn back to the case of Ad J\f = 2 QCD with massive hypermultiplets Then, there arise 
additional differentials corresponding to the derivatives of dS (2.28) with respect to masses. They are 
of the form ^ 

duj a = -- - _ - (2.33) 

It is again straightforward to use the residue formula^ 

dUjdLOjduK c j t jv' ! r . ■ , I o ! 

Fijk= res , {J, J,K, . . .} = {i,j,fc, ... a,P>7> ■••} (o "U1 



4 These punctures emerge as a degeneration of the handles of the hyperelliptic surface so that the A-cycles encircle 
the punctures. 

5 Let us note that, despite the differentials dw a have the pole at infinity, this does not contribute into the residue 
formula because of the quadratic pole of d\ in the denominator. 



and obtain 

J 1 v-^ Aijv c 



^ \ 1 6 v - 1 1 x -> 

^=Er + r + E x — ?E- 

_3_ | 2 | _J_ ^jN c ly 1 



^ ~ A^r + A, W + ^ * * " ' i¥=j ' 



A,Q ~ 2 U 

1 1 (2-35) 



1^1 1 1 

4 A Q/ 3 



These formulas immediately lead to the prepotential 

T = z E / (4) (««) - 1 E / (4) (^ - ™«) + i E z (4) (™« - ^) 



4^" A^ J v * " 16 



upon the identification A; = a,; and A Q = m a . 



(2.36) 



3 XJZ spin chain and 5d theories 

Now we are going to describe the generic integrable system behind the 5d theories - the inhomogeneous 
XXZ spin chain, and reproduce from the integrable system the main ingredients of the 5d SW theory: 
the spectral curve of and generating 1-differential dS. As in Ad case it allows to construct explicitly 
the holomorphic differentials on spectral curve and calculate the perturbative part of 5d prepotential. 

3.1 XXZ spin chain 

In 2 x 2 formalism the generalization of Ad formulas is straightforward: the Lax matrix for the XXZ 
(sl(2)) spin magnetic is defined on a cylinder and has the form |17| 

/ fie So - ur x e- So 2S- \ 

m = _ (3.1) 

y 25+ ne~ So - n~ x e So J 

In fact, from the point of view of integrable systems more natural spectral parameter is £ = log// 
so that the matrix elements of the Lax operator ([O]) are trigonometric functions of £. The Poisson 
bracket of the Lax operator (CO) 



{L(n) ? LQj>)} = [r(/x - J), m ® L(j/)] 



(3-2) 



is defined by the trigonometric r-matrix 



i 

r (£) = ~~ • i 7 ( CT i ® °"i + CT 2 ® f2 + cosh7r£er 3 ® cr 3 ) (3.3) 

sinn 7r£ 



so that (3.2) gives rise to the following Poisson algebra of Si's: 

{S+, S } = ±S±\ {S+, S-} = sinh2S (3.4) 
The second Casimir function is now 

K = cosh2S + 2S+S- (3.5) 



The non-linear relations (3.4) remind the commutation relations of quantum algebra U q (sl(2)), this is 



a generalization the fact that the XXX Poisson brackets (2.8) reproduce the classical algebra si (2). 



The value of deformation parameter q is inessential in (3.4) and can be put equal to unity; in fact, it 
is proportional to the radius of the target space circle R5 in corresponding 5d SUSY theory and can 
be easily restored replacing each generator Si by R$Si. Hereafter, we omit i?5 from all the formulas 
except for some special cases below. 

Again as in 4<i case one should consider the chain with N c sites, now with the Lax operators ( [j!l| ) 
associated to each site 

Li(fj,)%(ti) = % +1 (fi) (3.6) 

and commuting with each other for i ^ j; introduce the inhomogeneities which depend on the site 
of chain replacing £ — > £ — £j and impose periodic boundary conditions. As before the transfer matrix 
acting from i-th to i + iV c -tb. site is the product of the Lax matrices ( |2.5| ) 

T N M= ft Ih(£) (3.7) 

and the (quasi) periodic boundary conditions 

%+nM =-w*i(n) (3.8) 
provide the spectral curve equation now for the XX Z spin chain 

det(T Nc (n) +w ■ 1) = (3.9) 

The coefficients in this equation can be taken as a complete set of integrals of motion. 

The manifest form of the curve in the XXZ case can be derived using explicit expression for the 
Lax matrix fl3.1p and read^ 

w + 

2N, 

e 



' 2Af -(f) =2e-* r « g -I>P J v (e 2 Z)=2e- N ^-^(e 2N ^ + ... + e 2 J:^) (3.10) 

^+ 2 E 4 ^™ 1 V J 



Changing the variables e 2 ^ = fi 2 = A, w — * e S< ^ui, this curve can be rewritten in the hyperel- 
liptic form in A variables: 

w + 9lEM = 2 p NAXh Y* = Pl(\)-Q 2Nc (\), Y=Uv,-Q™M) (3.11) 
w 2 \ w I 



6 In fact, the coefficient in front of Pjv c i instead of 2, is equal to 2coshL| iVc e So i . It is, however, the integral of 
motion that can be put equal to any number. We always fix } t . 5o,i to be zero. 



while in terms of "true" spectral parameter £ this curve looks considerably more tricky. However, one 
can work with variable A taking into account that it lives on a cylinder or sphere with two marked 
points. In the equation ( |3.1l| ), the polynomial 



12N, 



(A) = J] (A 2 - 2 Kl e^X + e*) ^ J[ (a e^ +) ) (a - e*». H ) (3 . 12) 



defines the masses mj^ of the matter hypermultiplets. They turn out to be functions of the second 
Casimirs Ki of the algebra (3.4) at z-th site and corresponding inhomogeneity to be compared with 



(Ul) 



& ± -log ( Ki 



K'r 



(3.13) 



At the same time, the second polynomial Pn c (A) depends on the integrals of motion (mutually Poisson- 
commuting with each other) playing the role of moduli of gauge theory in the framework of the SW 
theory. Modulo these subtleties, the curve ( fr.ll ) is very similar to the curves arising for Ad theories. 
However, the difference coming from different generating 1-differentials dS turns to be very crucial 
(see, for example, the discussion of residue formula and perturbative prepotentials). 



It is easy to compare the curves (3.11) with those of |22| . rewriting them in the form 

IL sinh (£ - 



w - 



2 ]J sinh (f-Oi) 



(3.14) 



where we rescaled 



^ + ^2i^w, denoted the roots of the polynomial P Nc (\) ( |3.10|) A, = 



,2 a i 



and made use of formula (3.13) and the manifest form of the leading and the constant terms in this 
polynomial. Comparing (3.14) with ( |3.10|) , one finds, in particular, that 



6 



(3.15) 



The condition (3.15) is not, however, absolutely necessary in the context of XX Z chains and bd 
theories (we return to this point later in s.3.3). It emerges only in the standard XX Z chain. In 
the context of SW theory one needs rather to consider the twisted XXZ model |2l], [l^] []. It is 
characterized by the Lax operator 



fie 



s ,i-ti 



25- 



fie 



2S- ti 

-So.i-C« _ p^- 



(3.16) 



with generally non-unit constants o^, ft. These constants provide arbitrary coefficient in front of 



product in the r.h.s. of (3.14) (which can be also regulated by the integral of motion J^f So,i) and 



thus, break the condition (3.15). They are also important for careful treatment of dimensional trans- 
mutation (i.e. decreasing the number of massive hypermultiplets) Jl5| . The maximally degenerated 
case is related to the pure gauge theory when all the masses become infinite. In this case, one gets 
the relativistic Toda chain, first discussed in this context in || . 

There certainly exist all intermediate degenerations when all but Nf masses are brought to infinity. 
This system, corresponding to the gauge theory with Nf < 2N C massive hypermultiplets, is described 
by the chain with some of the sites degenerated and some of them not. The corresponding spectral 
curve has the form 



Qn } (A) 



N, 



= 2Pat c (A), Qjv,(A)=JJ(A 



N c 



^v c (A)=II( A - 



(3.17) 



7 We are grateful to S.Kharchev for the discussion of the twisted spin chains. 



or 

w + n^sinh(e-m a ) = 2e{Nc - Nfm+ J2 ai - 1/2 J2 ma ff sinh(c _ j (3 . 18) 
w J -- L 

i 

(where we again rescaled w — > (^/S+Zj m< *)iu; this curve also coincides with that in [p2[). Below we 
mainly discuss the generic degenerated system with Nf not necessarily equal to 2N C . As we already 
noted, the particular coefficient ai ~ 1 / 2 X m< * can be ,in principle, replaced by any other number, 
say put equal to unity, which literally corresponds to |2^| . It does not influence the result and we will 
discuss the reason below. 



Note that the form (3.14) of the spectral curve is perfectly designed for taking the 4c? limit. Indeed, 



one can restore /^-dependence in this formula multiplying each £j and mass parameter by R$. In 



terms of algebra (B.4) it means that each generator has to be multiplied by R5 and so does the spectral 



parameter. Then, one immediately reproduces the results of sect. 4. 
3.2 Relativistic integrable systems and 5d theories 

The general scheme of the SW anzatz presented above in the Ad case can be almost literally transferred 
to the "relativistic" hd M = 2 SUSY gauge models with one compactified dimension. One can 
understand the reason of this "relativization" in the following way. Considering four plus one compact 
dimensional theory one should take into account the contribution of all Kaluza-Klein modes to each 
4-dimensional field. Roughly speaking it leads to the 1-loop contributions to the effective charge of 
the form 

Tij ~ ^ l°g masses ~ ^ log ( a ij : + J ~ logJJ {Rsa-ij + m) ~ logsinhi? 5 a y (3.19) 

masses ra ^ ' m 

i.e. coming from 4c? to bd one should make a substitution — > smhR^aij, at least, in the formulas 
for the perturbative prepotential. The similar change of variables corresponds to relativization of 
integrable systems, which implies a sort of "Lie group generalization" of the "ordinary" integrable 
systems related to Lie algebras. In terms of 2 x 2 Lax representation, the "relativization" corresponds 
to replacing by trigonometric r-matrices and L-operators the rational ones: in algebraic terms this 
corresponds to turning from Yangians to quantum affine algebras. In particular, it means that now it 
is natural to consider both parameters £ — ^logA and z = logW as co-ordinates on cylinder, i.e. both 
bare spectral curves become the cylinders. 



Similar arguments imply that, instead of differential (2.20) dS^ — XdlogW, one has now to 
consider the differential ||, [T(| 

dS^ = ilogA^ (3.20) 

^ W 2 B W 

so that, despite the similarities between the 5g! and 4c? spectral curves, the periods of dS^ are different 
from the periods of dS^ . 

The difference with the 4c? case results also in a slightly different residue formula. Indeed, repeating 
the derivation of the paper , one can see that 

d^T „ dujjdwjdujK 

= 2 Fl -^rk^ (3.21) 



daidajdaji ~ fe dz^ 



The crucial distinction is that the contour in (3.21) encircles the "marked point" A = together with 



the zeroes of dz. Using this formula, we derive in next section the perturbative prepotential. 



Thus, one can interpret the 5d theory as the Ad theory with infinite number of (Kaluza-Klein) 
vector multiplets with masses M m = and infinite number of analogous (Kaluza-Klein) fundamental 
multiplets. Then, the presented 5c? picture naively corresponds to the infinitely covered "4c?" (see the 
spectral curve equation in ^-variable (3.14) living on a strip) also with infinitely many punctures. 
This, however, can be effectively encapsulated into hyperelliptic Riemann surface (3.11) with finite 
number of punctures, but, as a memory of infinitely many multiplets, the spectral parameter A now 
lives on a cylinder (on a sphere with two punctures). Meanwhile, the differential dS^ should be now 
of the form (3.20) which "remembers" its 4c? origin. 



Note that, in particular, the residue formula ( 3.21 ) can be written equally well in terms of variable 
£, instead of A. Then, in the denominator of (3.21) there arises the expression dzdt; and the residue 
is taken only at zeroes of dz (note also that this change of variables eliminates the factor 2 in ( 3.21 )). 
This form of the residue formula looks considerably more symmetric and we use its analog for the 6c? 
case in sect. 4. 



3.3 Perturbative limit of the SW anzatz 



In the perturbative limit, the hyperelliptic curve (3.11) turns into the rational curve 



w = 2 PnM = o nf c (A-A.) 



(3.22) 



with Xi = e 2ai , X)f c a i = 0; A Q = e 2m ° . We consider here the case of arbitrary number Nf > 2N C of 
massive hypermultiplets. First, let us construct the set of holomorphic differentials c?wj. 

The perturbative differential dS is 



dS = hogX dlog f P ^ 
2 \VW) 



or, equivalently, 



dS = £d£ [Y, coth (? - a i) ~ o coth (£ - m «) + Nc ~ o N f 



(3.23) 



(3.24) 



a=l 



The variation of ( 3.25 ) gives 



2SdS = -\og\d 



SX, 
X- Xi 



1 SXq \ 

2 ^ A - X a J 

a—l J 



dX SX 

X f-f Ai \X-Xt X-X Nc 

2—1 



1 c?A ^ SX a I X a 

2T^-XT U- K 

a—l x 



where it is implied that only jV c — 1 parameters Xi are independent so that , 
the set of the iV c — 1 independent canonical holomorphic differentials is 



(3.25) 



dX f X 
duji = —— 



X 



X\X-Xi A- A 



.v. 



XiN a dX 



(X-Xi)(X-X Nc ) 



i = l,...,N c -l, X, s X, X, (3.26) 



and one can check that the A-periods of g?5* are logAj = 2cii, while the condition of vanishing the 
sum of all residues gives the relation ( 3.15| ) J2i a i — \ m a- I n s.3.1, we discussed how this strange 



constraint should be removed off by using the twisted XX Z spin chain. The analytical explanation of 
this is that one can easily reproduce any lack of residues by putting some additional poles at infinities 
in £ or at A = 0, oo. It is analogous to the procedure in Ad. We shall see, however, that the same 
procedure does not work in d = 6. 



The additional differentials corresponding to the derivatives of dS (3.23) with respect to masses 
have the form |^| 

du n = -™ , " a . (3.29) 



dX X a 
2A A A a 



3.4 Perturbative prepotentials 



Now, using the manifest expressions for the differentials, we can calculate the perturbative prepoten- 
tials using the residue formulas. Let us first consider the case of pure SU(N C ) 5d gauge theory, when 
the prepotential is defined by the following residue formulas JIo| 

doJidujjdujk 



J~ ijk 2 



res 



<flogp=o dlogP(A)^ 



(3.30) 



since A = does not contribute into the residues. These residues can be easily calculated if one 
deforms the contour around the poles of du>iS. The final results have the following form 

6 COth CLiN c 



Tm = 2J cothaifc 



E 



T, 



"J 



cothajj + 4 coth + 2cothajAr c + 



coth a k N c , 
coth a k N c - N c , i^j, 



k^i,J,N c 

Tijk = 2 cotha;jy c - } cothai Nc - N c , i^j^k 

l^N c ljti,j,k,N c 

From these expressions it follows that the prepotential is given by formula 



(3.31) 



~Li 3 (e- 



h3 



3 

3 -ij 



;Li 3 (e 



— a-iajak 

i>j>k 



(3.32) 



Now let us turn to the computation of the perturbative part of prepotential for the 5d SU(N C ) theory 
with Nf massive fundamental hypermultiplets and compare the results with the answers appeared in 
literature (23]. In order to include massive hypermultiplets, we use the residue formula 

duidujjdujK 

(3.33) 



JK 



A = 



P(X) dX 1 
y/QW X 



{I, J.K, . . .} — {i, j, k, . . . | a, P, 7, . . .} 



°It can be instructive to describe hi detail why (2.21)) arises not only in Id, but also in 5rf models. Strictly speaking, 
in the 5d case, one should consider differentials on annulus - on sphere with 2 marked points. Therefore, instead of 
. d \ , one rather needs to take 



E 



coth (£ — ai)d^ - 



Taking now, instead of (2.29), the differentials (3.27) 

duji 



A + \i dX A + Ajv e dX 



A + Xi dX 
X- Xi 2X 



Aijv c dX 



A — Ai 2A A-Aiv c 2A (A — A;)(A — Ajv c ) 



(3.27) 



(3.28) 



we obtain again the formula (2.29). On the other hand, the mass differentials ( j3,29[ ) look different from (3.27). It turns 
out, however, that the difference does not contribute into the results and, therefore, one is free to choose any of these 
two possible mass differentials. We will return to this point in next subsection. 



and obtain 



r.Fm — E 



Ai + A 



AaT c ^ 1 Aj^Aq, 



A', 



= 1 + 2 



y 

AiAT c A,JV C ,~-f T ^kN, 



k^N, 
Aat c XiXjN, 



AfcAr 2 Aj Q Aiv a 

c a— 1 c 



E 

l—i,j,k 

A c 



AiJV c Ay 

^E 



2 h Nc 



2 , Ajv c a 
a— 1 c 



A a 



Aia AAr eC( 



A a 

V/V c a 



(3.34) 



« 7^ j 



Aia AAf c Q 



^ A 1 f \ 



2.A aa/3 



A, 



A a /3 



i.e. 



— E cotn a%k + 6cothaiAr c + cotha/yv c - 



k^N c 



+— y] (coth(aAr ( , — m Q ) — coth(ai — m Q )) 

a 

jFj i:) - = — cotha,y + 4cothaijv c + 2cothajjv c + y^ cothafejv c + - y^ coth(aAr c — m a )- 



N f 



J^ijfc = 2 E cotha;Ar c - E coth a iA r c + ^ E cotn ( a Ar c -m a )- 



l^i,j,k,N c 



(3.35) 



2-7 7 iia = coth(aj — TOq,) + coth(ajv c — m a ) — 2 

2^ r ija = coth(aAf ( , — m Q ) — 1 z 7^ j 
Zfiaa = (coth(a7v ( , — m a ) — coth(a i — ra a )) 

1T aaa = E coth ( a i ~ m a) + ^ E cot H m a - m t} ) - N c + ^J- + ^ 
4^00/3 = coth(m / 3 — m a ) — 1 

The "strange" numbers -1, -2, | that arise in the derivatives of the prepotential w.r.t. the mass moduli 
are due to the choice of the mass differentials in the form (3.29). Certainly, these numbers contribute 



1 o 



into at most quadratic part of the prepotential in gauge moduli. This does not effect the "physical" 
implications but does break the WDVV equations (l(J. Meanwhile, with the second possible choice 
(see footnote 5) 

dXX + X a d£ , . . , 

duj a = - — — = -— coth(£ - m a ) (3.36) 

these strange number are no longer present in the derivatives - and this turns out to be the only 
influence of the replace of the differentials. 

The difference between the two forms of mass differentials is due to different possible choices of 



the factor in the r.h.s. of (3.18). How it was already mentioned, this choice does not influence the 
result. It does, however, the WDVV equations. 

Thus, the bd perturbative prepotential can be written in the following plausible form 

T = i E /(5) m - 1 E /(5) ^ - w «) + Ye E / (5) k« - 

i,j i.a a,/3 

(3.37) 



where, as in (2.36) we use the functions f^ d \x) equal to 



/< 5 >(*) = ^ 3 -±Li 3 (e- 2 *) (3.38) 

4 XYZ chain and 6d theories 
4.1 General comments 

Now in the framework of general scheme developed in sect. 2 and applied to the 5d case in sect. 3, we 
are going to turn to the 6d case, or to the theory with two extra compactified dimensions, of radii R§ 
and Rq. The same argument as we used for the transition from four to five dimensions, i.e. taking 
account of the Kaluza-Klein modes, allows one to predict the perturbative form of the charges in the 
6d case as well. Namely, one should expect them to have the form 



(4.1) 



Ev-^ / m n 

log masses - ^ log I ay + — + — 

masses m,n ^ ^ ■ 

~ log JJ \ Rf>aij + m + n-jjrj ~ ^og6 \ R 5 a,ij i^J 

i.e. coming from Ad (5d) to 6d one should replace the rational (trigonometric) expressions by the 
elliptic ones, at least, in the formulas for the perturbative prepotential, the (imaginary part of) 
modular parameter being identified with the ratio of the compactification radii R$/R 6 . 

A similar change of variables would give rise to " elliptization" of an integrable systems. In terms 
of 2 x 2 Lax representation, the "elliptization" corresponds to replacing by elliptic r-matrices and 



-L-operators the rational (2.1) and trigonometric (3.1) one. In particular, it means that now it is 
natural to consider parameters £ = ^logA as a co-ordinate on torus, i.e. one of the bare spectral 
curves becomes a torus, while the other one remains to be a cylinder. 

In fact, we know the only system providing proper elliptization of the XXX spin chain, that is, the 



1 A 



XYZ Sklyanin chain G3| []. We propose it here as a candidate for the integrable system behind the 
6d theory. Below we describe the XYZ spin chain and demonstrate that it leads to the perturbative 
prepotential which coincides with our naive expectations (4.1), giving some support to our conjecture. 



4.2 XYZ chain 

Now we briefly describe the XYZ chain following || . All necessary facts on the elliptic functions can 
be found in Appendix. The Lax matrix is now defined on elliptic curve E(t) and is explicitly given 
by (see |l7j] and references therein): 

3 

L sk \o = s°i+t^Y, w ^) sa<j - ( 4 - 2 ) 

a=l 

where 

To keep the similarities with ||, we redefine in this section the spectral parameter £ — ► ijK, where 

K = — £ = f^n(O), k 2 ee so that if -> f as r -> ioo. This factor results into 

additional multiplier 7r in the trigonometric functions in the limiting cases below. 



The Lax operator ( |4.2| ) satisfies the Poisson relation (2.2) with the numerical elliptic r-matrix 
r(£) = i& J2a=i W a {£,) a a <S> <J a , which implies that S°, S a form the (classical) Sklyanin algebra [O: 



{S a ,S } =2i(?-) 2 (e b -e c )S b S c 

W (4.4) 



{S a ,S b } = 2iS°S c 

with the obvious notation: abc is the triple 123 or its cyclic permutations. 

The coupling constant — can be eliminated by simultaneous rescaling of the ^-variables and the 
symplectic form: 

S a = -S a S° = S° { , } -> -2^{ , } (4 5) 

g u 



Then 



L Skl (0 = S°l + i ]T W a (OS a a a (4.6) 

a=l 

= -i (e c - eb) S b S' 



{s a , 5°} 



a=l 

b 

(4.7) 



-iS u S 



c c 



One can distinguish three interesting limits of the Sklyanin algebra: the rational, trigonometric and 
double-scaling limits H. We are interested in only trigonometric limit here, since it describes the 



degeneration to the 5d case. In this limit, r — > +ioo and the Sklyanin algebra (4.7) transforms to 

{5 3 ,5°} = 0, {S\S } = iS 2 S 3 , {S 2 ,S } = -2iS 3 S 1 

~ ~ ~ ~ ~ ~ ~ ~ ~ ~ (4-8) 

{S\S 2 } = -iS°S 3 , {S 1 ,S 3 } = iS°S 2 , {S 2 ,S 3 } = -iS°S 1 

9 In [^J it has been already proposed that the XYZ spin chain should have something to do with the Nf = 2N C SW 
theory. We will see below, that this is indeed the case: the condition Nr = 2N C which can be easily broken in Ad and 
5d situations is very strict here. One may think of this theory as of blowing up all possible compactified dimensions in 
the Nf = 2N C theory with vanishing /3-function, see the discussion below. 



The corresponding Lax matrix is 

Lxxz = <5°1 

and r-matrix 



^— (sVi + S 2 a 2 + cosh^S 3 CT 3 



r (£) = — : — r (ci <8> cri + ct 2 ® 02 + cosh7r£(73 (g) cr 3 ) 
smn 7r£ 



(4.9) 



(4.10) 



Now one can note that the algebra ( |4.8| ) admits the identification So = e « +e~ h o , S3 — e~ — 
e s o 9 . With this identification and up to normalization of the Lax operator (4J]), we finally get the 
Lax operator (3^1) of the XX Z chain with the Poisson brackets (3.4). 

The determinant det2x2-^ Skl (£) i s equal to 

3 3 

det 2x2 i SkI (0 = So + J2 - = K ~ = K ~ M * x ( 4J1 ) 

a— 1 a— 1 

where 

K = S 2 Q +J2e a (r)S 2 a M 2 = £^ (4.12) 

a— 1 a— 1 

are the Casimir operators of the Sklyanin algebra (i.e. Poisson commuting with all the generators S , 
S\ S 2 , S 3 ). 



Now, in order to construct the spectral curve, note that the determinant of the monodromy matrix 



T N M = T\N c > l >i LSk ^-^ 



is 



q(o = det 2x2 T Nr (o = n det 2x 2 i ski (e = n ( Ki - M ?p(t - ^)) 



n^l g '(^- m c) n TT 1 

- — — 7^-7 yo — 11 A 



1=1 



(4.13) 



nr=i^(e-^) 2 



where we used formulas ( A.8), ( A. 9 ) and ( A.llj) from Appendix in order to rewrite K{ — M 2 p(£ — £j) 



9.(f- m | + >)8.(g-m|-)) 

as -/i?' ^. j • riere 



At 



Mf[9'Mf (±) 
6i{Ki) ' m * 



£2 i /^i 



and Ki is a solution to the equation 

6l{Ki) KiOKO) n 2 0i(0)[6t(0) + 6j(0)) 



02 («i) M 2 [^(0)f 



3[^.(0)r 



Again 



(4.14) 



(4.15) 



(4.16) 



As a particular example, let us consider the case of the homogeneous chain || (all = in ( |4.13| )) 
TrTjv c (£) is a combination of the polynomials: 

P(0 = Pol^ [x) + (s), (4.17) 

where [^M is integral part of and the coefficients of Pol^ and Pol^ are Hamiltonians of the 
XFZ model. The spectral equation (3.£) for the XYZ model is: 



= 2P(0 



(4.18) 



1 / ■ 



where for the homogeneous chain P and Q are polynomials in x = p(£) and y — |p'(£)- Eq. ( 4.1S| ) 
describes the double covering of the elliptic curve E(t): with generic point £ £ E(t) one associates the 



two points of T, XYZ , labeled by two roots w± of equation ( 4.18 ). The ramification points correspond 
to w+ = w_ = ±V<2, or Y = | (w - = ^P 2 - Q = 0. Note that, for the curve ( ff.lSD , a; = oo is 
noi a branch point, therefore, the number of cuts on the both copies of E(t) is N c and the genus of 
the spectral curve is N c + 

In the case of a generic inhomogeneous chain the trace of monodromy matrix looks more compli- 
cated, being expression of the form 



Tr [] f 

iV c >i>l 



Ski 



{a,} 2V >i>l 



(4.21) 



Af c >i>l 



a 4 = 0, aj 



1,2,3; W (0 = 1 



It is clear from ( |4.3| ) that each term in the sum (4.21) (and thus the whole sum) has simple poles in 
all inhomogeneities £ = Thus the trace of the monodromy matrix can be represented as 



(4.22) 



53 a * = X) & 



■>z_. "■" (4-23) 
where cij should be identified with (shifted) gauge moduli - see discussion in sect. 3.1. The equation 



( 4.18 ) can be finally rewritten as 



UaiiO*(Z~rn a ) 

*n£°i0.(e-&) a 



Qo 



n&f.«-fi) 



(4.24) 



We see that in 6d case it is necessary always to require Nf = 2N C in order to have a single-valued 



meromorphic functions on torus. Again we also meet the constraint (4.23). In the 5d case, it was 
removed, along with the condition Nf — 2N C twisting the spin chain. Since no twisted XYZ chain is 
known we are forced to preserve both constraints, their origin will be discussed below q 



The curve (4.24) has no longer that simple expression in terms of polynomials of x and y similar 



to (4.17). Thus, in the inhomogeneous case it is not convenient to work in terms of "elliptic" variables 
x and y (analogous to A in the 5d case) instead we choose to deal with the variable £ <E C/I\ 



10 Rewriting analytically H XYZ as a system of equations 

3 



Y[(x - e„), 



the set of holomorphic 1-diffcrcntials on £ 

dx x a dx 
v = — , V a = 

y yY 



y = l l ( > - • i- y 

a=l 

XYZ can be chosen as 



Q 



o, 



c^dx 
Y 



(3 = 0,. 



(4.19) 



(4.20) 



with the total number of holomorphic 1-differentials 1 + ( [^M + l) + ( [ ^"2 3 ] + ■*■) = N c + l being equal to the genus 
of T, XYZ . 



11 Formula (4.23) implies that, while the sum of all the would-be gauge moduli a; still remains a constant, it is no 



longer zero. Therefore, one would rather associate with the gauge moduli the quantities a; shifted by the constant 
^2 m a . We do not know the proper field theory interpretation of this condition. 



Given the spectral curve and integrable system one can immediately write down the generating 
1-diffcrcntial dS. It can be chosen in the following way 

d3 XYZ s _ logu< (4.25) 
w 

Return to our oversimplified example of the homogeneous chain. Now, under the variation of moduli 
(which are all contained in P, while Q is moduli independent), 

S(dS XYZ ) - ^ = - Sm = %5P (4.26) 



and, according to (4.20), the r.h.s. is a holomorphic 1-diffcrcntial on the spectral curve ( 4.1q ). 

The residue formula in the 6c? case is of the form completely analogous to the 4d case (and to the 
5d case, when working in terms of the variable £) 

d 3 T duJidujjduJK tA ~x 

res : — — (4-27) 



daidajdax dz=o dzd£ 
and we use it in the next subsections to derive the perturbative prepotential 



4.3 Perturbative limit of the XYZ spin chain 

A generic XYZ curve ( i . 24| ) in the perturbative limit turns into 



(4.28) 



where ho = 2 f° . Then, the rescaled generating differential ( 4.25 ) is 

V Qa 

dS = [J2 log^* (£ - Oi) - \J2 l °Z * ^ - "*»)) « + logftod* = 
^C((-a j )--^CK-m 1 



(4.29) 



where the linear piece disappears due to ^ a,; = i ^ m Q . This expression should be compared with 



( p.24| ) in the 5d case. 



Note that the sum of residues vanishes only provided the constraint ( 4.23 ) is fulfilled. It can not 
be removed by adding an additional pole at infinity, since the variable £, in contrast to the Ad and bd 
cases, lives on the compact surface with out boundaries (torus) and, therefore, is always finite. 

Now one can take the total variation of dS (we neglect the trivial variation Sho giving rise to the 
holomorphic differential d£ on the bare torus) 



N c j N f 

SdS XYZ £* gc(£ - aJScndt - - C(£ - m a )5m a <% 



(4.30) 



i=l 



where C(£l r ) i s defined in Appendix (see (A6)). Therefore, the differentials related to gauge moduli 



arc 



dun = (C(? - Oi) - C(£ - a Nc )) d£ 



(4.31) 



while those related to masses are 



1 



duj a = - - (C(£ - m a ) ~ C(C ~ ojvJ) d£ 



(4.32) 



The second term in this expression is due to the condition (4.23). Thus, one can see that, in the 6d 
case, mass moduli are practically identical to the gauge ones. 

The residue formula acquires the form 

dujidtojdiOK 



Tijk = res dl =0 



7^ " u d^d\og^= 



(4.33) 



4.4 Perturbative prepotentials 



Having the explicit formulas for the differentials ( 4.31 ) and ( 4.32 ), one can use the residue formula 
( 4.33| ) to calculate the perturbative prepotential in the 6d case. The result is 



Fin = ^2((aik) + 6C(ffliiv c ) + ^2 ^( afc ^)" 

k^i k^N c 



a 

Fuj = + 4 C(aiAfJ + 2((a jNc ) + ^2 C( a kN c ) + - ^2C(aN c - m a ), i ^ j 



k^i,j,N a 



T ijk = 2 C( 



aiN c 



C( a iN c ) + ^^2C{aN c - m a ), i^j^k 



l^N c l=jti,j,k,N a 

2Fiia = C( a i -m a ) + ((a Nc -m a ) + C( a N e i) + ^2 Ct"*^) ~ \ ~£( aN <= ~ m ^ 



JV, 



N f 
1 «rA-, 



— — — v— ^— 1 . ■ - 

2^ a = C(a Nc - m a ) + ((a Nci ) + ((a Ncj ) + ^ C( a N c k) — ^ ^ ajv <= ~ m ^ i ^ J 

fc p 

l( Nc - l Nf - \ 

\k p J 

1 - 1 ( Nc - 1 Nf - \ 

i /3^a \ fc /3 / 



2-7 r i CtQ = ((cti - m a ) + -C(oiJV c ) , } 



(4.34) 



3 7v \ 

+ ^C(w Q - a Nc ) 



1-. 



1 / ^ - 1 Nf - 

,) + 5 EW- E?( 



JV C _ 1 W/ _ 

4^ Q/ 3 = C(™a - a» c ) + C( m p - ojvJ + CKjvJ - X! £( a JV<=fc) ~ 2 ^2 ^ aN - ~ m ^ 



k 



%F a p~t = -C( m a -aN c )- C(mp -a Nc )- CKjvJ + ^ C(aN a k) ~ ^ ^2 C( a N c ~ m p ) 



These formulas imply that the perturbative prepotential is given by 

T = i E / (6) ko - z E / (6) ( a * - m «) + E / (6) ( to « - ^) 



v v/ 4^ J v ' 16 

i,j i,a ot,p 



(4.35) 



where a, i = 5 Zi m a an< ^ the function (x) defined so that 

f^"{x)=\o g 9 if {x) (4.36) 
can be expressed through the elliptic tri-logarithm Li3. g (a;) |2J]: 



(4.37) 



3 4 120 2 



5 Discussion 

Now we are going to discuss the results. We have reproduced the spectral curves for the SUSY SU(N C ) 
gauge theories with Nf massive hypermultiplets in d — 4, 5 and 6 dimensions. Using these curves and 
the corresponding differentials, we have calculated the perturbative parts of the prepotentials in all 
these cases and the results have the same form. 

Results. The spectral curve in all cases can be written in the form 

W + Q^H& = 2 pW(o (5-1) 

w 

or 

™ ■ 1 - 2P(d) ^ , W = ^4= (5.2) 
The generating differential dS is always of the form 

dS = £d\ogW (5.3) 
The perturbative part of the prepotential is always of the form 

T = I E - I E f (d) ^ ~ m «) + Yq E - (5.4) 

i,j i,a a. /3 

The concrete forms of the functions introduced here are: 

No N c N c . . 

i i j *VS SV 

(in P^(£)i there is also some exponential of £ unless 2V/ = 27V C , see ( |3.18[ )) 

/«>(*) = Ags, /<«(*) = yj/W (* + = i |*»| - ±Li 3 (e- 2 l"l) , 

/ m «=£/ w (*+i+f)=E/ <5) Ht)= M 



l , 

3 



x 3 l 



-Li3 )9 2 ' a: '^ + quadratic terms 



on 



so that 

fW" = ]ogp, f&"(x) = logsinhx, fM"(x)=]og0.(x) (5.8) 

Note that, in the 6d case, Nf is always equal to 2N C . The variables & above are inhomogcneities in 
integrable system, and, in d — 5,6, there is a restriction = J~)$j — \ X) m a which implies that 
the gauge moduli would be rather associated with et^ shifted by the constant j^- m a . 



On perturbative limit. Now let us make some comments concerning the perturbative limit in all 
these cases. In fact, the procedure is quite clear in the Ad and 5d cases: this limit corresponds to 
big values of gauge moduli (and masses) compare to Aqcd and can be taken as -> 0, which is 



exactly the perturbative limit in asymptotically free theories. Looking at (5.2), one can see that this 
procedure corresponds to £ going to infinity (as ^qq D ) and W going to infinity or to zero depending 
on the choice of a particular list of hyperelliptic spectral curve From the point of view of the 
effective charges TJy it separates logarithmic terms from the power (instantonic) contributions. 

This reasoning certainly fails when d — 6. Indeed, since £ G C/T it can not be taken to infinity. 
In terms of prepotential it means that one can not distinguish between \og9 terms and the other since 
^-function is restricted. However, there is a way out of this problem, although a little bit artificial. 
Namely, instead of moving £ one may just turn to zero the coefficient in front of the product in Q^(£). 
The perturbative contribution to the 6d prepotential calculated in the paper correspond exactly to 
this limit. 

As it is well-known, 5d and 6d theories are non-renormalizable. Therefore, the status of the 
perturbative contributions has to be understood better (see also [^3[ ||, |i|). We are not going 
to discuss this question in detail here, instead let us point out that the problems of renormalizability 
are related only with maximally quadratic part of the prepotential getting stronger than logarithmic 
divergencies in 5d and 6d perturbative theory. Our approach does not control these terms and the 
properties of (compactified) bd and 6rf theories wc discuss do not depend on quadratic contributions. 

The first property is the periodicity of the bd prepotential under — > + 2mnR^,~ 1 , with integer 
n, which is evidently correct up to quadratic terms (being changed due to the shifts in cubic terms). 
This shift is induced by changing the sheet of the logarithm in dS^ ( 3.20D (either logA, or logw, 



depending on the chosen form of dS^) as one can easily obtain from (2.25). It explains, in particular 



why the ambiguity in definition of dS^ 5 " 1 ( 3.20| ) does not influence the result. Similarly, in 6d case 



one needs double-periodicity to guarantee the unambiguous prepotential despite the multi-valuedness 



of generating differential dS^ ( 4.25 ) . This is indeed the case: /( 6 ) in ( 4.37 ) is a double-periodic 
function. 

The second property of the prepotentials is typically to satisfy the WDVV equations |)| |l(| . This is 
strongly related to integrable structure of the theory and should have a topological origin. The WDVV 
equations are equations for the third derivatives and do not depend on quadratic terms. On the other 
hand, these equations requires to consider masses and gauge moduli on equal footing. Therefore, the 
WDVV equations are sensible to the constant corrections to the third derivatives of the prepotential, 
even if some of these derivatives are taken w.r.t. masses. 

This issue has a lot of to do with the proper definition of mass differentials (see sect. 3), we are 



12 In 4d W — » oo as Agp^ " . However, one can renormalize the gauge moduli, masses, W and £ so that they remain 

finite as Aqcd 0. It results into the coefficient Aqq^, ' in front of the term ^ in (5.2). The analogous procedure 
can be also performed for the other sheet, when W — > 0. In 5d the dependence on K.qq]j can not be reduced to the 
only coefficient, this is an indication of non-renormalizability of 5d theory. 



going to return to this problem elsewhere. 



5d and 6d theories — cubic terms, Nf — 2N C etc. Next point we are going to discuss is the cubic 
(in moduli) terms appearing only in higher-dimensional theories. First, take the 5d case, where these 
terms corresponds to the Chern-Simons term in the field theory Lagrangian Tr(^4 A F A F). In our 
computation they appear due to the constant piece (N c — \Nf) in the brackets in (3.24). 



In fact, let us choose the special values of the second Casimir functions at m sites so that the 
polynomial Q(X) acquires the form \ m Q Nf (X), i.e. m of 2N C factors in (3.12) turn into A m . In this 



case one gets the constant piece N c -m-Nf/2 in dS {5 \ with an integer m < 2N C — Nf. Therefore, 



the coefficient in front of cubic term in ( 3.37 ) can be made equal to h(N c — m — Nf/2) 



Restoring the dependence on i?5 in ( |3.37 ) (at — > aiR$ and m a — > m a R^) one can study the different 



limits of the system. The simplest limit corresponds to Ad case and is given by i?5 — > 0. In this limit, 



f( 5 \x) — > f^(x), the cubic terms vanish and we reproduce the perturbative 4d prepotential ( 2.36 ) 



At the level of integrable system it is enough to replace Si — > RsSi, [i — > e R& ^ in the Lax operator 



( p.l| ) ( (|3.16| )) in order to reproduce the Lax operator of the XXX (higher sl(j>)) spin chain, see |lfj 
for details). 

Another interesting limit is the limit of flat 5d space-time, i.e. R§ — > oo. In this limit, only cubic 



terms survive in the prepotential ( 3.37| ) (one should carefully fix the branch of f^(x) which leading 
to appearing of the absolute value in ( |5.9| )): 

_ 1 

■ mp\ + 

(5.9) 




In fact, there are two different sources of cubic terms p^] . The first one is from the function f^(x). 
Since this function can be obtained as the sum of the 4d perturbative contributions over the Kaluza- 
Klein modes, these cubic terms have perturbative origin and come from the 1-loop (due to the well- 
known effect of generation of the CS terms in odd-dimensional gauge theories). 

The second source of the cubic terms is due to the bare CS Lagrangian. As it was shown in |p3| , 
one can consider these terms with some coefficient c c ;: 

Ccl \ ^ 3 

lL a > (5.10) 

i 

restricted only to satisfy the quantization condition c c i + -J- G Z an d the inequality |c c /| > N c J-. 



In our formulas, we easily reproduce these conditions so that (5.9) coincides with p3[, provided the 



prepotential is defined in a fixed Weyl chamber. We leave, however, the wall crossing jumps out of 



the discussion. Note that the curve ( 3.11 ) with this arbitrary m coincides with that proposed in J22|. 

Note also that it is clear why the WDVV equations are not necessarily satisfied when m ^ 0. 
Since it corresponds to degeneration of an integrable system, one fixes some moduli (masses) - the 
parameters to vary in the WDVV equations (see jl0| for details). 

In the 6d case the calculations are quite similar, and one would expect the similar cubic terms 
correspond to Tr(F A F A F). However, despite all the difference in calculations is due to slightly 



different mass differentials (3.36) and (4.32) the bare cubic terms in 6d are absent. This occurs due 
to Nf — 2N C (since the coefficient in front of cubic terms is proportional to N c — hNf) and, therefore, 



cancelation of the constant term in the brackets of ( 4.29| ) as compared with ( |3.24 ). This is one of the 
problems with 6c? theories: one can not deform (twist) the theory so that the constraints Nf = 2N C 
and ( 4.23 ) look unavoidable. We do not understand clearly their meaning at the level of field theory. 
As for the "quantum generated" terms, i.e. those coming from the function they are certainly 
presented, with the coefficient equal to that in the bd case. 



Finally, let us point out that the Lax matrices (2T), (3J) and (4J2) (and analogously though 
more involved for the transfer matrix (4.21)) can be possibly determined from a linear differential 
equations DLi(£) ~ #( 2 )(£ — £j) (a Id Hitchin). These equations could be interpreted as arising in 
the framework of D-brane interpretation of SUSY gauge theories [|l2j [HJ . It might indicate that the 
standard Hitchin-like systems and the spin-chain like models, in spite of distinctions in the formulation 
of the Poisson structures etc |fj"7f , can be treated within some unified approach in the framework of 
higher-dimensional SUSY gauge theories and corresponding brane configurations of string theory. We 
are going to return to this problem elsewhere. 



Appendix 

The definition of elliptic functions we use throughout the paper slightly differs from standard defini- 
tions that can be found, for example, in p7| . In this Appendix we define all the functions we need 
and write down manifestly some useful identities necessary for sect. 4. 

We define the Weierstrass p-function to be the sum 

-Uqq -(-OO 

e(*l T )= E 7TZ \ v ~ E , 1 ^ =72+g(g 2 ) (A- 1 ) 

(f + m + nT) z ^— ' (m -I- nrY f z 

m,n— — oo v ' m,n— — oo v ' 

where means omitted term with m = n = 0. Thus defined p-function is double periodic, with 
periods 1 and t = ^ and differs from the standard definition by the factor (2lj)~ 2 and rescaling 
£ — > 2ui£. According to this definition, the values of p(£|r) in the half-periods, e a = e a (r), a = 1,2, 3, 
also differ by the factor (2a;) -2 from the standard definition. 

The complex torus E(t) can be defined as (D /Z © r Z with a flat co-ordinate £ defined modulo 
(1,t). Alternatively, any torus (with a marked point) can be described as the elliptic curve 

y 2 = (x-e 1 )(x-e 2 )(x-e 3 ), x = p(0, V = \p\£.) (A.2) 



and the canonical holomorphic 1-differential is 

at = 2— (A.3) 

y 

Other elliptic functions used in the paper are ^-function 

C(Ck) = 7+ E it — - — + —^ + ( r ^ 2 )-7 + ^ 3 ) (a.4) 

F — ' V F — m — nr m + nr (m~\-nr) z / F 

^ m,n=— oo x v / * 



and ^-functions 



+00 



0i(£|T) = M£|T) = <9,(eir) = i E (-!) r 



i(n-£) 7tt «7r(2n-l)£ 



n— — 00 



e 2 (e|r) = e 01 K|r)= 2 el 



(n-±y-KT e i-K(2n-l)t 



+ OO 



(A.5) 



n— — 00 



+00 



n— — 00 



There exist the following relations connecting different elliptic functions and constants with 1 
and 6>-constants 0(0) = 0(0\t): 

where 



f , 10*' (0) 



with 



p(0 = -C'(0 = e 3 + 
es(r) = 



y.(o)^(0 



3 L 2 



'4(0) 0,(0 

(O) + 0f(O)] 



To simplify some formulas, we also introduce the logarithmic derivative of the 0-function 

C(£|T)=C(£|T)-7?(TX = [log0*(O]' 

In sect. 4, we also use the identity for 0-functions |^7j 

9,{x + y)9*(x - y)el(Q) = 9l{x)6l(y) - 9j(x)9t(y) 



-functions 
(A.6) 

(A.7) 
(A.8) 

(A.9) 
(A.10) 
(All) 
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